Pb-Pb collisions at y/s^N = 2.76 TeV in a multiphase transport model 
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The multiplicity and elliptic flow ol charged particles produced in Pb-Pb collisions at center 
of mass energy ^/snn = 2.76 TeV from the Large Hadron Collider are studied in a multiphase 
transport (AMPT) model. With the standard parameters in the HIJING model, which is used as 
initial conditions for subsequent partonic and hadronic scatterings in the AMPT model, the resulting 
multiplicity of final charged particles at mid-pseudorapidity is consistent with the experimental data 
measured by the ALICE Collaboration. This value is, however, increased by about 25% if the final- 
state partonic and hadronic scatterings are turned off. Because of final-state scatterings, particular 
those among partons, the final elliptic flow of charged hadrons is also consistent with the ALICE data 
if a smaller but more isotropic parton scattering cross section than previously used in the AMPT 
model for describing the charged hadron elliptic flow in heavy ion collisions at the Relativistic 
Heavy Ion Collider is used. The resulting transverse momentum spectra of charged particles as well 
as the centrality dependence of their multiplicity density and the elliptic flow are also in reasonable 
agreement with the ALICE data. Furthermore, the multiplicities, transverse momentum spectra 
and elliptic flows of identified hadrons such as protons, kaons and pions are predicted. 

PACS numbers: 25.75.-q, 12.38. Mh, 24.10.Lx, 24.85. +p 
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Experimental data on Pb-Pb collisions at ^/sjvn = 
2.76 TeV from the Large Hadron Collider (LHC) by the 
ALICE Collaboration have recently become available [J- 
13]. For most central collision bins (0 — 5%), the mid- 
pseudorapidity density of charged particles was found to 
be 1584 ± 4(stat.) ± 76(sys.) Qj, which is a factor of 2.2 
higher than that observed in central Au-Au collisions at 
v / s]vlv = 200 GcV from the Relativistic Heavy Ion Col- 
lider (RHIC). This multiplicity density was reproduced 
by the HIJING2.0 model with a more modern set of par- 
ton distribution functions [f| and has helped to constrain 
the gluon shadowing parameter in the model Q. Fur- 
thermore, the elliptic flow in non-central collisions was 
found to have values similar to those in collisions at 
RHIC energies Q. According to Ref. Q, the similarity 
between the elliptic flows at LHC and RHIC is consis- 
tent with the predictions of the viscous hydrodynamic 
model. In Ref. Q, this similarity has further led to the 
conclusion that the specific viscosity rj s /s of the quark- 
gluon plasma produced in heavy ion collisions at LHC 
has a similar value as in heavy ion collisions at RHIC. 
The observed elliptic flow at LHC can also be described 
by a kinetic model with a short-lived parton stage but 
strong final-state hadronic scatterings [9(. As shown in 
Ref. [lCj . adding final-state scatterings, which are essen- 
tial for generating the observed elliptic flow [ll|, to the 
HIJING model as implemented in a multiphase transport 



(AMPT) model [12[ would reduce the predicted charged 
particle multiplicity at mid-rapidity as a result of their 
appreciable contribution to the longitudinal work [l3j . It 
is thus of interest to study both the multiplicity and el- 
liptic flow in Pb-Pb collisions at ^Jsnn = 2.76 TeV by 
using the AMPT model. 

This paper is organized as follows. In Sec. [HI wc briefly 
review the AMPT model and discuss its parameters. We 
then study in Sec. Illll thc multiplicity and in Sec. HVI thc 
elliptic flow in Pb-Pb collisions at ^/snn = 2.76 TeV. 
Finally, we give some discussions in Sec. [V] and the con- 
clusions in Sec. IVII 



II. THE AMPT MODEL 

The AMPT model is a hybrid model with the ini- 
tial particle distributions generated by the HIJING 
model [HI of the version 1.383. In the version of string 
melting, which is used in the present study, hadrons pro- 
duced from the HIJING model are converted to their va- 
lence quarks and antiquarks, and their evolution in time 
and space is then modeled by the ZPC parton cascade 
model [ll| with the differential scattering cross section 
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In the above, t is the standard Mandclstam variable for 
four momentum transfer, a s is the strong coupling con- 
stant, and (i is the screening mass in the partonic mat- 
ter. After stopping scattering, quarks and antiquarks 
are converted via a spatial coalescence model to hadrons, 
which arc followed by hadronic scatterings via the ART 
model [ll. 
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In previous studies of heavy ion collisions at RHIC, 
some of the parameters in the HIJING model were var- 
ied in order to reproduce the measured charged particle 
multiplicity. Specifically, instead of the values a = 0.5 
and b = 0.9 GeV~ 2 used in the HIJING model for 
the Lund string fragmentation function f(z) oc — 
z)° exp(— b m^/z), where z is the light-cone momentum 
fraction of the produced hadron of transverse mass m± 
with respect to that of the fragmenting string, the val- 
ues a = 2.2 and b = 0.5 GcV~ 2 were used in the AMPT 
model [Tol | • Also, the values a s = 0.47 and /i = 1.8 
or 2.3 fm , corresponding to a total parton scattering 
cross section of 10 or 6 mb, were used for the parton 
scattering cross section in the AMPT model to describe 
measured elliptic flow [T3, E3 and two-pion correlation 
functions [19j . In the present study of heavy ion colli- 
sions at LHC, we again vary the values of these parame- 
ters to fit measured data. As shown below, a reasonable 
description of the measured charged particle multiplicity 
density at mid-pseudorapidity is achieved if the default 
HIJING values of a = 0.5 and b = 0.9 GcV~ 2 arc used in 
the Lund string fragmentation function. Also, a smaller 
QCD coupling constant a s = 0.33 and a larger screen- 
ing mass [i — 3.2 fm" 1 , corresponding to a smaller (1.5 
mb) but more isotropic parton scattering cross section, 
are needed to reproduce the observed elliptic flow. 



III. RAPIDITY DISTRIBUTIONS AND 
TRANSVERSE MOMENTUM SPECTRA 



model with string melting using different values for the 
Lund string fragmentation parameters and parton scat- 
tering cross sections. It is seen that the values of a = 2.2 
and b = 0.5 GeV -2 , that correspond to a larger string 
tension in the Lund string fragmentation function, gives 
a larger multiplicity density at mid-pseudorapidity than 
that from the default HIJING values of a = 0.5 and 
b = 0.9 GcV~ 2 that correspond a smaller string tension 
in the fragmentation function. The multiplicity density 
is, however, not sensitive to the parton scattering cross 
section a. For the smaller value of tr = 1.5 mb, which is 
required to reproduce the measured elliptic flow at LHC 
as shown later in Sec. IIV1 the resulting charged parti- 
cle multiplicity of 1,536 ± 4 at mid-pseudorapidity ob- 
tained with the default HIJING parameters for the string 
fragmentation function is consistent with the measured 
data by the ALICE Collaboration. In the following, we 
will thus use, unless stated otherwise, those parameters 
that correspond to the smaller values of string tension 
and parton cross sections in the AMPT model for our 
study. We note that in obtaining the relation between 
the centrality c and the impact parameter b, we have 
used the empirical formula c = 7rb 2 /a- ln [20l | with the 
nucleus-nucleus total inelastic cross section a- m ~ 784 
fm 2 calculated from the Glauber model. Also, we have 
used in our analysis about 1 , 000 AMPT events for the 
multiplicity density, transverse momentum spectra, and 
the total elliptic flow, and about 10,000 events for the 
transverse momentum dependence of the elliptic flow. 
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FIG. 1: (Color online) Pseudorapidity distribution of charged 
particles in Pb-Pb collisions at ^/sjvjv = 2.76 TeV and central- 
ity of — 5% from the AMPT model with string melting for 
larger string tension (filled pentagons) and for smaller string 
tension with a parton scattering cross section of 10 mb (filled 
stars) or 1.5 mb (filled squares). 

We first show in Fig. [T] the pseudorapidity distribution 
of charged particles in the most central 5% (b < 3.5 fm) 
Pb-Pb collisions at ^/IWn = 2.76 TeV from the AMPT 
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FIG. 2: (Color online) Pseudorapidity distribution of charged 
particles in Pb-Pb collisions at ^/snn = 2.76 TeV and cen- 
trality of - 5% from the AMPT model with both final-state 
partonic and hadronic scatterings (filled squares), with only 
partonic scatterings (filled triangles), and without final-state 
interactions (filled circles). 

We have also studied the effect of final-state interac- 
tions (FSI) on the charged particle multiplicity. As shown 
by filled circles in Fig. [21 turning off the final-state par- 
tonic and hadronic scatterings in the AMPT model en- 
hances the charged particle pseudorapidity distribution. 
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The charged particle multiplicity at mid-pseudorapidity 
is now 1925 ± 5, which is essentially the value from 
the HIJING1.0 model Q and is about 25% higher than 
that with the inclusion of final-state scatterings, which is 
shown by filled squares. The value is reduced to 1642 ± 4 
after including partonic scatterings, and this shows that 
the charged particle pscudorapidity distribution is much 
more sensitive to partonic than hadronic scatterings. 
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FIG. 3: (Color online) Dependence of charged particle pseu- 
dorapidity density per half participant (dN c h/ drj) / ({N part } /2) 
on the number of participants in Pb-Pb collisions at ^/snn = 
2.76 TeV from the AMPT model with string melting and from 
the ALICE data @|. 



Figure |3] displays the dependence of the charged parti- 
cle multiplicity density at mid-pseudorapidity < 0.5) 
per 1/2 participant, (dN c h / drj) / ((A part ) /2), on the num- 
ber of participants. It is seen that results from the AMPT 
model with string melting (open circles) are roughly con- 
sistent with the experimental data, although for periph- 
eral collisions the values are higher than the ALICE data. 

For identified hadrons such as protons, kaons, and pi- 
ons as well as their particles and antiparticles, their ra- 
pidity distributions in the — 5% centrality of same col- 
lisions arc shown in Fig. 01 The multiplicity ratio of mid- 
rapidity protons, kaons and pions is roughly 1 : 2 : 16. 
We note that the multiplicities of each species are again 
larger if final-state scatterings are not included. 

In Fig. [5j the transverse momentum (pr) spectrum 
of mid-pseudorapidity (|?/| < 0.8) charged particles (left 
panel) and those of protons, kaons and pions (right panel) 
are shown. It is seen that the AMPT model describes rea- 
sonably the charged particle px spectrum from the AL- 
ICE data for small pr but gives smaller values for large 
Pt- For identified hadrons, the abundance of pions is the 
largest at all pt except at intermediate pr where the pro- 
ton yield becomes comparable. A similar phenomenon 
was observed in heavy ion collisions at RHIC 



21|, and 



it was attributed to the enhanced production of baryons 
from quark coalescence f22T - |24j ] . 
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FIG. 4: (Color online) Rapidity distributions of protons, 
kaons and pions in Pb-Pb collisions at ^/snn = 2.76 TeV and 
— 5% centrality from the AMPT model with string melting. 
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FIG. 5: (Color online) Transverse momentum spectrum of 
mid-pseudorapidity (|?j| < 0.8) charged particles (left panel) 
as well as those of protons, kaons, and pions (right panel) in 
Pb-Pb collisions at ^/snn = 2.76 TeV and — 5% centrality 
from the AMPT model with string melting. Experimental 
data from the ALICE Collaboration Q| are shown by solid 
squares in the left panel. 



IV. ELLIPTIC FLOWS 

In the present study, we determine the ellip tic flow 
using the two-particle cumulant method [25], [If 



v 2 {2} = V(cos(2A0)), 



(2) 



where Ac6 is the azimuthal angular difference between 
particle pairs within the same event and (• • • ) means av- 
erage over all possible pairs. The error of the elliptic flow 
is calculated by using the method in Ref. [26| . 

In Fig. [SI we show by filled circles the transverse mo- 
mentum dependence of the elliptic flow «2{2} of mid- 
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FIG. 6: (Color online) Transverse momentum dependence 
of the elliptic flow obtained from the two-particle cumulant 
method for charged particles in Pb-Pb collisions at ^/snn = 
2.76 TeV and 40 - 50% centrality from the AMPT model 
with string melting using a parton scattering cross section of 
1.5 mb (filled circles) or 10 mb (filled stars). Corresponding 
experimental data from Ref. @] are shown by filled squares. 
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FIG. 7: (Color online) Centrality dependence of the elliptic 
flow obtained from the two-particle cumulant method for all 
charged particles (upper panel) and for protons, kaons and 
pions (lower panel) in Pb-Pb collisions at ^/snn = 2.76 TeV 
from the AMPT model with string melting. The ALICE data 
(filled squares) are from Ref. @|. 



DISCUSSIONS 



pseudorapidity (|?/| < 0.8) charged particles at central- 
ity of 40 - 50% (10.0 fm < b < 11.2 fm) from the 
AMPT model and compare it with the ALICE data (filled 
squares) Q. It is seen that the elliptic flow from the 
AMPT model is consistent with the experimental data 
at pt < 2.5 GeV/c but is larger at p T > 2.5 GeV/c. The 
latter is likely due to an overestimated non-flow effect in 
the AMPT model. We note that a satisfactory under- 
standing of elliptic flow of high- p-^ p articles in heavy ion 
collisions at RHIC is still lacking [27ll28j . Using the larger 
parton scattering cross section of 10 mb, corresponding 
to a larger strong coupling constant (a s = 0.47) and a 
smaller screening mass (/i = 1.8 fm -1 ) as used in de- 
scribing the elliptic flow at RHIC, would lead to a larger 
elliptic flow (filled stars) than the ALICE data. 

The centrality dependence of the i>2{2} of mid- 
pseudorapidity charged particles with transverse mo- 
menta 0.2 GeV/c < pt < 5 GeV/c is shown in the upper 
panel of Fig. [7] by using again the smaller parton scat- 
tering cross section of 1.5 mb. It is seen that the results 
from the AMPT model (open circles) are consistent with 
the experimental data (filled squares) except for central 
and peripheral collisions where they arc slightly larger 
and somewhat smaller, respectively. The smaller elliptic 
flow for centralities larger than 40 — 50% is likely related 
to the softer pt spectrum in the AMPT model. In the 
lower panel of Fig. [71 we compare the centrality depen- 
dence of the elliptic flows of protons, kaons and pions. 
Interestingly, the elliptic flow is larger for kaons than for 
pions. 



The above results from the AMPT model for heavy 
ion collisions at LHC were obtained with different val- 
ues for the parameters in the Lund string fragmentation 
function and in the parton scattering cross section from 
those used for heavy ion collisions at RHIC. The values 
of a = 2.2 and b = 0.5 GeV" 2 used at RHIC for the Lund 
string fragmentation function correspond to a string ten- 
sion that is about 7% larger than that for a = 0.5 and 
b = 0.9 GeV -2 used in the present study. A smaller 
string tension at LHC is consistent with the fact that the 
matter formed at LHC is hotter and denser than that 
at RHIC. The higher temperature reached at LHC also 
makes it possible to understand, according to the lattice 
results [23,[3(3, the smaller QCD coupling constant and 
the larger screening mass needed to describe the elliptic 
flow at LHC. 

We can compare the property of the quark-gluon 
plasma produced in heavy ion collisions at LHC with that 
at RHIC by considering the ratio of its shear viscosity and 
entropy density, i.e., the specific viscosity. In the kinetic 
theory, the shear viscosity is given by -q s = 4(p} / (15a tr ) , 
where (p) is the mean momentum of partons and Ctr is 
the parton transport or viscosity cross section defined by 



[ dt—(l - cos 2 
J dt y 



187ro^ 



1 
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1 + y 2 /E 2 
jj?/E 2 



- 2 



,(3) 



where E is the center of mass energy of the colliding 
parton pair. In obtaining the second line of above equa- 
tion, we have used Eq. ([T]). For quark-gluon plasma of 
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massless quarks and gluons at temperature T, we have 
(p) = 3T and E ~ y/l8T. The entropy density s of 
the quark-gluon plasma, which is modeled by quarks and 
antiquarks of current masses in the AMPT, is simply 
s = (e + P)/T = 4e/(3T) = 96T 3 /tt 2 if only up and down 
quarks are considered. The resulting specific viscosity is 
thus 



3tt 



1 



40a 2 



(9 + ^)m(^-)-i8 



(4) 



The initial temperature of heavy ion collisions in the 
AMPT model can be estimated from the average energy 
density of mid-rapidity partons at their average forma- 
tion time, which is about 46.0 GeV/fm 3 at LHC and 19.5 
GcV/fm 3 at RHIC. Using e = 72T 4 /tt 2 for the baryon- 
free quark and antiquark matter, we obtain an initial 
temperature of about 468 MeV at LHC and about 378 
MeV at RHIC. The value of r] s /s in the early stage of 
the quark-gluon plasma formed in these collisions, when 
most of the elliptic flow is generated, is thus about 0.273 
at LHC and 0.085 (for /i = 1.8 fm" 1 ) or 0.114 (for fx = 2.3 
far 1 ) at RHIC. 

We note that for fixed values of /z and a s , which is 
the case in the AMPT model, the value of rj s /s as given 
by Eq. ((4]) increases as the temperature of the quark- 
gluon plasma decreases. In a more realistic description, 
the screening mass depends on temperature, i.e., /i = 
gT with g 



\J Akol s [31|. In this case, the ratio r\ s js 
becomes independent of temperature. It is then of great 
interest to extend the AMPT model to determine the 
screening mass from the QCD coupling constant and the 
local temperature of the partonic matter as in Ref. (32j 
and use the resulting model to study if using the same 
QCD coupling constant a s or same specific viscosity r] s /s 
in the AMPT model would describe the elliptic flow in 
heavy ion collisions at both RHIC and LHC, as in the 
findings of Rcfs. 0,Q based on the hydrodynamic model. 



y/sNN — 2.76 TcV. We have found that the measured 
multiplicity density and elliptic flow of charged particles 
at mid-pseudorapidity can be reasonably described by 
the model if the parameters in the Lund fragmentation 
function are taken to be those used in the default HIJING 
model and that a smaller but more isotropic parton scat- 
tering cross section than that used for heavy ion collisions 
at RHIC is used. As at RHIC, the final-state partonic and 
hadronic scatterings were found to be important as they 
would reduce the charged particle multiplicity density at 
mid-pseudorapidity by about 25%. The smaller parton 
cross section needed to describe the measured elliptic flow 
at LHC than at RHIC has led to a larger estimated 7? s /s 
in the quark-gluon plasma produced at LHC than that at 
RHIC. However, this result may be due to the use of con- 
stant screening mass in calculating the parton scattering 
cross section. Taking into account the medium depen- 
dence of the screening mass may reduce this difference. 
Furthermore, the transverse momentum spectra and the 
centrality dependence of the multiplicity and elliptic flow 
have been studied and they are roughly consistent with 
the ALICE data. We have also made predictions for the 
multiplicity distributions and elliptic flows of identified 
hadrons such as protons, kaons, and pions. Comparisons 
of these predictions with future experimental data will 
further enhance our understanding of heavy ion collision 
dynamics at LHC and the properties of produced quark- 
gluon plasma. 
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